Nodal gap structure in Fe-based superconductors 
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To understand the origin of the nodal gap structure realized in BaFe2(As,P)2, we study the 
three-dimensional gap structure based on the three-dimensional ten-orbital Hubbard model with 
quadrupole interaction. In this model, strong spin and orbital fluctuations develop by using the 
random-phase-approximation. By solving the Eliashberg gap equation, we obtain the fully-gapped 
s-wave state with (without) sign reversal between hole-like and electron-like Fermi surfaces due 
to strong spin (orbital) fluctuations, so called the s±-wave (s ++ -wave) state. When both spin 
and orbital fluctuations strongly develop, which will be realized near the orthorhombic phase, we 
obtain the nodal s-wave state in the crossover region between s++-wave and s±-wave states. The 
obtained nodal s-wave state possesses the loop-shape nodes on electron-like Fermi surfaces, due to 
the competition between attractive and repulsive interactions in fc-space. In contrast, the SC gaps on 
the hole- like Fermi surfaces are fully-gapped due to orbital fluctuations. The present study explains 
the main characters of the anisotropic gap structure in BaFe2(As,P)2 observed experimentally. 

PACS numbers: 74.20.-z, 74.20.Fg, 74.20.Rp 



I. INTRODUCTION 

Since the discovery of Fe-based high-T c superconduc- 
tors by Kamihara et al.[l], their many-body electronic 
properties have been studied very intensively. Figure 1 
shows a typical phase diagram of Fe-based superconduc- 
tors. In the under-doped regime, the second-order or- 
thorhombic (O) structure transition occurs at Ts, and 
the stripe- type magnetic order is realized at Tn < T$. 
In the O phase, the orbital polarization n xz ^ n yz is re- 
alized, where n xz r yz \ is the filling of d xz ^ yz ^ orbital [2]. 
Also, sizable softening of shear modulus C§§ [3-5] and 
the renormalization of phonon velocity [6] indicate the 
development of orbital fluctuations near the orthorhom- 
bic phase. Strong spin fluctuations are also observed near 
the magnetic ordered phase. 

These main characters of the phase diagram should 
be understood in order to clarify the mechanism of su- 
perconductivity. Theoretical studies of orbital polariza- 
tion had been proposed in Refs. [7-9]. However, non- 
magnetic structure transition cannot be explained based 
on the Hubbard model, once we apply the mean-field 
approximation or random-phase-approximation (RPA). 
To solve this problem, we had recently improved the 
RPA by including the vertex correction (VC) for the 
susceptibility that is dropped in the RPA [10]: By ap- 
plying this self-consistent VC (SC-VC) method to the 
Hubbard model for Fe-based superconductors, both spin 
and orbital fluctuations mutually develop, and both the 
O structure transition and the softening of Cqq can be 
explained. Note that the "electronic nematic state" with 
large in-plane anisotropy of resistivity or magnetization 
well above Ts [11-13] also indicates the occurrence of the 
(local) orbital order [14]. 

The phase diagram in Fig. 1 indicates that both spin 
and orbital fluctuations could be closely related to the 
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FIG. 1: (Color online) A typical phase diagram for Fe-based 
superconductors. Ts is the structure transition temperature, 
which is expected to be induced by orbital polarization n xz > 
n yz according to the ARPES measurements and the sizable 
softening of Cee- Tjv is the magnetic transition temperature. 



mechanism of high-T c . Up to now, the spin-fluctuation- 
mediated s±-wave state [15-18] and orbital-fluctuation- 
mediated s ++ -wave state [19, 20] had been studied based 
on the multiorbital models. The s±-wave state has the 
sign reversal of the gap between hole-like Fermi surfaces 
(h-FSs) and electron-like Fermi surfaces (e-FSs) , whereas 
the s_|_+-wave state has no sign reversal. Experimentally, 
the robustness of T c against impurities in many Fe-based 
superconductors [21-24] indicates the realization of the 
s ++ -wave state, at least in "dirty" compounds with high 
residual resistivity [25]. Also, the "resonance-like" hump 
structure in the neutron inelastic scattering [26] can be 
explained if the s ++ -wave state is realized. 

Although fully-gapped s-wave state is realized in many 
optimally-doped high-T c compounds, nodal s-wave state 
(accidental node is not protected by symmetry) is also 
observed in some compounds with lower T c [27]. The 
appearance of the accidental node strongly indicates the 
presence of "competing pairing interactions" [28] . In the 



2 



spin fluctuation scenario, if the xy-orbital hole-pocket 
is under the Fermi level, spin fluctuations in xz + yz- 
orbitals develop at Q = (tt, 0) while those in xy-orbital 
develop at Q' = (tt,tt/2), and this frustration gives a 
nodal gap structure around the sy-orbital part on the 
e-FS [29]. This mechanism was studied in detail in Ref. 
[30] by introducing phenomenological pairing interaction. 
However, the xy-orbital h-FS presents in real compounds. 
In this case, fully-gapped s±-wave state is obtained by 
the RPA since spin fluctuations develop at Q — (tt, 0) in 
all d-orbitals, consistently with neutron experiments [31]. 

Very interestingly, in optimally-doped BaFe2(As,P)2, 
nodal gap structure with high-T c (~ 30K) is realized. 
The superconducting (SC) gaps on h-FSs are fully- 
gapped and almost orbital-independent [32, 33], consis- 
tently with the orbital fluctuation scenario in Ref. [20]. 
Also, loop-shape nodes on the e-FSs are observed by 
angle-resolved thermal conductivity measurement in the 
vortex state [34] and ARPES measurements [32, 33]. Al- 
though different gap structure was reported in Ref. [35], 
small density-of-states around the loop-node is consistent 
with the small Volovik effect in the specific heat measure- 
ment [36, 37] and other bulk measurements [38]. These 
results indicate the existence of competing pairing inter- 
actions, and the study of these facts would be significant 
to understand the mechanism of high-T c superconductiv- 
ity. 

In this paper, we theoretically study the origin of 
the nodal gap structure in BaFe 2 (As,P) 2 , in order to 
obtain a significant information of the pairing mecha- 
nism of Fe-based superconductors. For this purpose, we 
construct the three-dimensional (3D) ten-orbital tight- 
binding model for BaFe2(As,P)2, and calculate the dy- 
namical spin and orbital susceptibilities due to the com- 
bination of Coulomb and quadrupole interactions. By 
solving the Eliashberg gap equation, we obtain the fully- 
gapped s±-wave (s++-wave) state due to strong spin (or- 
bital) fluctuations. When both spin and orbital fluctua- 
tions strongly develop, which will be realized near the O 
phase, nodal s-wave state with loop-shape nodes on the e- 
FSs is realized due to the competition between attractive 
and repulsive interactions. It is realized during a smooth 
crossover between s++- and s±-wave states [19, 39]. In 
contrast, the SC gaps on the h-FSs are fully-gapped due 
to orbital fluctuations. Thus, the present study explains 
the main characters of the gap structure in BaFe2(As,P)2- 

In Refs. [19, 20, 40-42] present authors have shown 
that small quadrupole interaction induced by Fe-ion os- 
cillations gives rise to the large antiferro- and ferro- 
orbital fluctuations. In addition, we had developed the 
spin+orbital fluctuation theory in multiorbital Hubbard 
model by including the VCs to the susceptibilities, which 
are neglected in the RPA [10]. It was found that the 
Azlamasov-Larkin type VC due to Coulomb interaction 
produces large effective quadrupole interaction. The 
emergence of the orbital fluctuations due to the VC 
is also recognized in a simple two-orbital model, using 
the self-consistent VC method [43] as well as newly de- 



veloped two-dimensional renormalization group method 
(RG+cRPA method) [44]. 
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FIG. 2: The Fermi surfaces in the (a)fc z = plane and (b)k y = 
plane of the present ten orbital model for the filling n = 6.0. 
The green, red, blue and black lines correspond to xz, yz, xy 
and z 2 orbitals, respectively. In (b), there are three h-FSs 
(FSl, FS2 and FS3) and four e-FSs (FS4 and FS5). 



II. FORMULATION 

In this paper, we set x and y axes parallel to the near- 
est Fe-Fe bonds, and the orbital z 2 , xz, yz, xy, and 
x 2 — y 2 are denoted as 1, 2, 3, 4, and 5 respectively. First, 
we perform the local-density-approximation (LDA) band 
calculation for BaFe 2 As 2 and BaFe 2 P 2 using WIEN2K 
code based on the experienced crystal structure. Next, 
we derive the ten-orbital tight-binding model that repro- 
duces the LDA band structure and its orbital character 
using WANNIER90 code and WIEN2WANNIER inter- 
face. [45] Using the obtained two sets of tight-binding 
parameters (hopping integrals and on-site energies), the 
parameters of BaFe2(Asi_ a; P a ;) 2 are well approximated 
by making a linear combination of them with a ratio of 
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1 — x : x [46]. In this paper, we use the tight-binding 
parameters for x = 0.30. The obtained kinetic term is 
given as 

ttO _ ,aa,b/3 f 

n — l lm C laa,<T C mbf3 ,a 

abafilma 

= E ^c^-^'O^w w 

abafilmcj k 

where a, b represent the unit cell, a, (3 (= A, B) represent 
the two Fe sites, l,m= 1 — 5 represent the d orbital, and 
a = ±1 is the spin index. R a , a is the position of Fe- 
site, cj is the creation operator of the electron, and 

with a = 6 and a = (3 (a ^ b or a ^ j3) is the local 
potential (hopping integral). 

However, the rry-orbital h-FS given by the LDA is too 
small compared to the experimental results by ARPES 
measurements. In order to increase the size of xy-orbital 
h-FS, we introduce the following orbital-dependent po- 
tential term around the r-point: 
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where e; is the energy shift of the orbital-/ at r-point. We 
put e xy = 0.02 eV, e xz — e yz — —0.01 eV and the others 
are 0. The FSs in this model are composed of three h-FSs 
around r-point and four e-FSs around X- and F-points. 
Figure 2 show the obtained FSs in the (a) k z — and 
(b) k y — planes, respectively. The electron filling per 
Fe-site is n = 6.0. In Fig. 2 (b), there are three h-FSs 
(FS1, FS2 and FS3) and two c-FSs (FS4 and FS5). We 
call FS4 (FS5) the outer (inner) e-FS. 

Next, we explain the interaction term. We introduce 
both the Coulomb interaction (U, U', J = (U - U')/2) 
and quadrupole interaction. The latter are induced by 
the electron-phonon (e-ph) interaction due to Fe ion os- 
cillations as follows, [41] 
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where gi{ui) = gi ■ Wp/(w, 2 +w^), and g t = g^O) is the 
quadrupole interaction at co n = 0. cjd is the cutoff energy 
of the quadrupole interaction. Or is the quadrupole op- 
erator and V^ ua d has many non-zero off-diagonal elements 
as explained in Refs.[19] . As explained in Ref. [19], g\ 
(52) is induced by in-plane (out-of-plane) Fe-ion oscilla- 
tions. In this paper, we put g\ = gi = g unless other- 
wise noted. Also, the Azlamasov-Larkin type VC due to 
Coulomb interaction produces large effective quadrupole 
interaction g\ [10]. Thus, the quadrupole interaction in 
eq. (3) is derived from both the VC and e-ph interaction. 



Now, we perform the RPA for the present model, by 
using 32 x 32 x 16 k meshes. The irreducible susceptibility 
in the ten orbital model is given by 
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where q = (q,u>i), k = (k,e n ), e n = (2n + l)irT and 
uji = 2lnT are the fermion and boson Matsubara fre- 
quencies. G{k) = [ie n + fi — /ifc ln ] _1 is the d electron 
Green function in the orbital basis, where h^ n is the ma- 
trix elements of H kln and fi is the chemical potential. 
Then, the susceptibilities for spin and charge sectors in 
the RPA are given by [47] 
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J , l\ = I4 7^ l 2 = h 

0. otherwise 



where a,fl = A, B. 

In the RPA, the enhancement of the spin susceptibility 
X s is mainly caused by the intra-orbital Coulomb inter- 
action U, using the "intra-orbital nesting" of the FSs. 
On the other hand, the enhancement of \ c in the present 
model is caused by the quadrupole-quadrupole interac- 
tion in eq. (3), utilizing the "inter-orbital nesting" of the 
FSs. The magnetic (orbital) order is realized when the 
spin (charge) Stoner factor a s i c ), which is the maximum 

eigenvalue of T s ^x^°\q, 0), is unity. When n = 6.0, the 
critical value of U is U CI = 1.18 eV, and the critical value 
of g is (? cr = 0.23 eV for [7 = 0. Hereafter, we set the 
unit of energy as eV. 

Next, we explain the linearized Eliashberg equation. In 
order to obtain the fine momentum dependence of the SC 
gap, we concentrate on the gap functions only on the FSs 
as done in Ref. [29]: We used 40 x 16 k points for each 
Fermi surface sheet. In the presence of dilute impurities 
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("•imp <C 1), the linearized Eliashberg equation is given 

as [29] 

Z 4 (fc,e„)A E A*(fc,e„) 

FS 



^3 EE/ ^V«(Jb,fc',e B -e m ) 



7rT 
(27 



/FSj 
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where Ae is the eigenvalue that reaches unity at T = T c . 
i and j denote the FSs, and A*(fe, e„) is the gap function 
on the i-th FS (FSi) at the Fermi momentum k. The 
integral in eq. (10) means the surface integral on FSj. 
The paring interaction V in eq. (10) is 

V^(k,k',e n -e m ) = Ut iati (k)Uup,i(k) 



v = v c + v s + v {0 \ 

2 V ; 



(12) 
(13) 
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where Ui a ^(k) — (k;la\k;i) is the transformation uni- 
tary matrix between the band and the orbital represen- 
tations. 

In cq. (10), Z is given as 



Z\k,e n ) = 1 + 



(15) 



where is the impurity induced quasiparticle damp- 
ing rate. Here, we calculate the damping rate using T- 
matrix approximation. We consider the case of Fe-site 
substitution, where the impurity potential / is diagonal 
in the d-orbital basis. The T-matrix for an impurity at 
a(= A or B) site is given as 



T a (e n )= i-I a G? oc {e n ) 



(16) 



which is fc-indcpendent in the orbital basis. Here, If v — 

Idi.i> is the impurity potential, and G" oc is the local Green 
function given as 

[Gioc]^(e„) = ^E G «'( fe '' e «) 
fe' 

'-^Zl^f^W^'), (17) 
where s n = sgn(e„). 



In the T-matrix approximation, which is exact for 
"imp "C 1, the normal self-energy in the band diagonal 
basis is given as 

SK(k,e n ) = n im pE^*a,i(W'(cn)Ui'a,i(fc), (18) 

ll'a 

where n lmp is the impurity concentration ratio. Then, 
the quasiparticle damping rate is given as 



7 i (fc,e„) = -Im^Ej 1 (fe,e„)s r . 



(19) 



Also, SY^I is the impurity-induced anomalous self-energy 
given as 
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In this calculation, we simplify the energy dependence 
of V. We assume that (£ = c, s) can be separated into 
the momentum and orbital dependent part V^(k, u>i = 0) 
and energy dependent part g^(coi): 



Ve(k,w l ) = V*(k,wi = 0)xgt(u l ). 



(22) 



We calculated V^(k,coi = 0) without approximation. On 
the other hand, g^(uJi) is determined approximately as 



X4ax(w/ = o) 



(23) 



where V,| ax (0) is the largest value of Vf '£^ l4 (k,coi = 0) 
for any a,/3,li, and fe. It is verified that this simplifi- 
cation affects the momentum dependence of the SC gap 
functions only quantitatively, although the obtained Ae 
is quantitatively underestimated. Thus, this approxima- 
tion would be appropriate for the present purpose, that 
is, the analysis of the anisotropy of the SC gap. 



III. SUPERCONDUCTING GAP 

In this section, we analyze the linearized Eliashberg 
equation, eq.(10), using the 3D model of BaFe2(As,P)2 
for n = 6.0. Hereafter, we use 32 x 32 x 16fc meshes for 
calculating charge and spin susceptibilities. We assume 
that J = J' and U = U' + 2J, and fix the ratio J/U = 
1/6. In solving the Eliashberg equation, we used 40 x 16fc 
points for each Fermi surface sheet and 512 Matsubara 
frequencies. In this paper, we perform the calculation 
T = 0.005 and cj d = 0.02. 
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A. s±-wave SC gap mediated by Spin fluctuations 
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FIG. 3: (Color online) Obtained SC gap functions for U = 
1.15 and g — 0. (a),(b) SC gap functions on the h-FSs in 
k z — and k z = n planes. The green, red, blue and black 
lines correspond to xz, yz, xy and z 2 orbitals, respectively. 

(c) k z dependence of the SC gaps on the h-FSs in k y = 
planes. Horizontal node appears on the FS3 around k z = ±7T. 

(d) ,(e) SC gap functions on the e-FSs in k z = and k z = tt 
planes. (f),(g) 3D gap functions on the outer and inner e-FSs. 



First, we study the spin-fluctuation-mediated s±-wave 
superconducting state for U < U CI by putting 5 = and 
«im P = 0. Here, we put U = 1.15 (a s = 0.98), and 
the obtained eigenvalue is Ae = 1.01. The obtained gap 
structure is almost independent of a s . First, we discuss 
the SC gaps on the h-FSs. Figures 3(a) and (b) show the 
obtained gap functions on the h-FSs in the k z — and 
tt planes, respectively. The definitions of 9 and FS1-5 
are shown in Fig. 2. In the k z = plane, the SC gap 
size weakly depends on the orbital character of the FSs. 
However, in the k z — tt plane, the SC gap size strongly 
depends on the d-orbital. Especially, the SC gap on the 
z 2 -orbital FS is almost zero and negative, reflecting the 
small spin fluctuations in the z 2 -orbital because of the 
absence of the intra z 2 -orbital nesting. (Note that z 2 - 
orbital is absent on the e-FSs.) The horizontal node is 
clearly recognized in the SC gap in the k y — plane 
shown in Fig. 3 (c). The obtained horizontal node on 
FS3 near k z — tt is consistent with the previous RPA 



calculation by Suzuki el al [46]. 

The obtained horizontal node would contradict to the 
four-fold symmetry of the thermal conductivity [34] and 
the small Volovik effect in the specific heat measure- 
ment [36, 37]. According to ARPES measurements, the 
horizontal-node was reported in Ref. [35] , whereas it was 
not observed in Refs. [32, 33]. 

Next, we discuss the SC gaps on the e-FSs. Figures 
3 (d) and (e) show the obtained gap functions on the 
e-FSs in the k z = and tt planes, respectively. As we 
can see, line nodes do not appear on the e-FSs. This 
result is consistent with the analysis in Ref. [31], that is, 
the s±-wave gap on the e-FSs is fully-gapped if the h-FS 
made of xy-orbital appears. Note that the SC gaps for 
k z = tt in Fig. 3 (e) are obtained by rotating the gaps in 
the k z = plane in (d) by tt/2. Also, Fig. 3 (f) and (g) 
show 3D gap functions on the outer and inner FSs (FS4 
and FS5), respectively. On both e-FSs, the SC gap on 
the "flat part" is larger than that on the "high curvature 
part" . 



s++-wave SC gap mediated by orbital 
fluctuations 



Next, we study the orbital-fluctuation-mediated S-|— (-- 
state superconducting state for g < g cr by putting U = 
and riimp = 0. Here, we put g — 0.22 (a c = 0.98), and 
the obtained eigenvalue is Ae = 0.59. The obtained gap 
structure is almost independent of a c . Figures 4 (a) and 
(b) show the obtained gaps on the h-FSs in the k z = 
and 7T planes, respectively. In highly contrast to the spin 
fluctuation scenario, the gap size on z 2 -orbital FS is com- 
parable with that on the other FSs, since strong orbital 
correlations are developed in all d-orbitals: Note that 
the quadrupole interaction possesses many non-zero in- 
terorbital matrix elements. The present numerical result 
is consistent with our previous calculation using the 2D 
5-orbital model. [20] 

Figure 4 (c) shows that the SC gap size of each h- 
FSs is approximately independent on k z , which is consis- 
tent with the small orbital dependence of the SC gap 
in (Ba,K)Fe2As2 and BaFe2(As,P)2 observed in Refs. 
[32, 33]. Figures 4(d) and (e) show the obtained gaps 
on the e-FSs in the k z = and tt planes, respectively. 
Figure 4 (f) and (g) show the 3D SC gap functions on 
the outer and inner e-FSs (FS4 and FS5), respectively. 
Thus, the obtained SC gaps on the e-FSs is isotropic for 
any k z . 

We also discuss the SC gap functions in the case of g\ — 
g and 172 = in eq. (3). Figure 5 shows k z dependence 
of the SC gaps on the h-FSs for g = 0.24 (a c = 0.98) and 
U = 0. In this case, the gap function on the z 2 -orbital 
h-FS is smaller compared to the case of g\ ~ <?2 = 9 in 
Fig. 4. On the other hand, the obtained SC gaps on 
the e-FSs are almost isotropic, similarly to the results of 
9i = 92 = 9- 
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FIG. 4: (Color online) Obtained SC gap functions for g = 0.22 
and (7 = 0. (a),(b) SC gap functions on the h-FSs in k z = 
and kz = it planes. The green, red, blue and black lines 
correspond to xz, yz, xy and z 2 orbitals, respectively, (c) k z 
dependence of the SC gaps on the h-FSs in k y = plane. 
Used colors are same as (a) and (b). (d),(e) SC gap functions 
on the e-FSs in k z — and k z = it planes. (f),(g) 3D gap 
functions on the outer and inner e-FSs. 

gap function 
on the h-FSs (ky=0) 




0.005 0.01 
A 

FIG. 5: (Color online) Obtained k z dependence of the SC 
gaps on the h-FSs in k y — plane. Used parameters are 
gi = 0.24, g 2 = and U = 0. 



C. Loop-shape node due to the competition of spin 
and orbital fluctuations 



Recently, several measurements observed the nodal gap 
structure in BaFe2(Asi_ a; P a: )2 [32-34]. This compound 



is very clean, and very accurate measurements of gap 
structure have been performed. They present a signifi- 
cant challenge for theories to reproduce the observed gap 
structure. However, as discussed in subsection III-A and 
B, we cannot reproduce the line- nodes on the electron 
FSs when either spin or orbital fluctuations solely de- 
velop. 



spin 
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k 7 =n 
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FIG. 6: (Color online) Formation of the nodal s-wave gap 
(shown in Fig. 7) due to the competition of orbital fluc- 
tuations (=inter-orbital attraction) and spin fluctuations 
(=intra-orbital repulsion). Green, red, blue lines correspond 
to xz, yz, xy, and z 2 -orbitals, respectively. 

Here, we study the emergence of highly anisotropic s- 
wave state due to the strong orbital and spin fluctuations. 
In the phase diagram of BaFe 2 (As 1 _ a; P : j;)2, both T/v and 
Ts decrease to zero at almost the same critical point 
x c R3 0.3. This fact means that both spin and orbital 
fluctuations become comparable in magnitude at x ~ x c . 
Here, we consider the case that the s++-wave state is re- 
alized by stronger orbital fluctuations. As increasing the 
spin fluctuation with momentum Q, A& and A^+q are 
suppressed when both k and k + Q are on the FSs with 
the same orbital character, and finally the sign change 
Afc-Afc + Q < could be achieved. Such strong anisotropy 
originates from the competition between the attractive 
interaction of V c and repulsive interaction of V s in eq. 
(13). As shown in Fig. 6, strong spin fluctuations on the 
xy-orbital (due to intra rry-orbital nesting) produce the 
loop-shape node on the e-FS. Similar "anisotropic s-wave 
gap modified by the spin fluctuations" is considered to 
be realized in (Y,Lu)Ni 2 B 2 C.[28] 

Hereafter, we present numerical results in the presence 
of small amount of impurities (1=1 and nj mp = 0.03), just 
to make the SC gap functions smoother. Figure 7 show 
the results of nearly s ++ -wave state with nodal structure 
on outer e-FS. We put g = 0.204 and U = 1.011, (a c = 
0.980, a s = 0.859), and the eigenvalue is Ae = 0.50. 
Figures 7 (a)-(c) show the obtained SC gaps on the h- 



FSs in the k z = plane, k z = n plane, and k y 







plane, respectively. The obtained SC gaps on the h-FSs 
are nearly isotropic and orbital- independent, similarly to 
the results in Fig. 4. Especially, the gap size of the z 2 - 
orbital h-FS is large even in the presence of loop-shape 
nodes on e-FSs. 

Figures 7 (d) and (e) show the obtained SC gaps on 
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FIG. 7: (Color online) Obtained SC gap functions for g — 



0.204, U = 1.011 and n- m 



0.03. (a),(b) SC gap functions 



on the h-FSs in k z = and k z — ty planes, (c) k z dependence 
of SC gaps on the hole FSs in k y — plane. (d),(e) SC gap 
functions on the e-FSs in k z — and k z — it planes. (f),(g) 
3D gap functions on the outer and inner e-FSs. The green 
lines represent the gap nodes. 



the e-FSs in the k z — and n planes, respectively. The 
SC gap on the inner e-FS is fully opened, and its sign 
is same as that on the h-FSs. On the outer e-FS, in 
contrast, the SC gap shows the sign change near # = 0, 7r 
(0 = 7r/2, 37r/2) in the k z — plane (k z — it plane). 
This sign change is caused by strong spin fluctuations in 
the xy-orbital, as we have explained in Fig. 6. In this 
case, the SC gaps on the h-FSs remains fully-gapped, due 
to the fact that the band-mass of h-FSs is larger than 
that of the e-FSs. As results, closed loop-shape nodes 
appear in the flat part on the outer e-FS, as recognized 
in Figs. 7 (f) and (g). This gap structure is consistent 
with the prediction given by the angle-resolved thermal 
conductivity under the magnetic field [34]. 

As increasing the value of U (or reducing n; mp ) slightly, 
the area of the sign reversed part on the outer e-FS in- 
creases, and the SC gap on the inner e-FS also shows 
the sign reversal. Figures 8 (a) and (b) show the SC 
gap functions on the e-FSs for g = 0.204 and U — 1.017 
(a c = 0.980 and a s — 0.864). The obtained eigenvalue 
is Ae = 0.50. The obtained gap functions are approxi- 
mately given by shifting the gaps in Figs. 7 (d) and (e) 
downwards, and line nodes appear on both the inner and 



inner e-FS 

\ outer e-FS 




(c)outer e-FS (FS4) 
^ k z =;c 

VflPf I ' 001 



0.001 



-0.001 



k 7 =-iz 



(d)irmer e-FS (FS5) 
k 7 =;t 



k 7 =-K 



FIG. 8: (Color online) Obtained SC gap functions for g = 
0.204 and U = 1.017. (a),(b) SC gap functions on the e-FSs 
in k z — and k z = n planes. (c),(d) 3D gap functions on 
the outer and inner e-FSs. The green lines represent the gap 
nodes. 



outer e-FSs. As described in Figs. 8 (c) and (d), closed 
nodal loops appear in the flat part on the outer e-FS and 
in the high curvature part on the inner e-FS. The SC 
gaps on the h-FSs are almost the same as those shown in 
Fig 7 (a)-(c), so we do not show them. 



inner e-FS 




(c)outer e-FS (FS4) 
k 7 =it 



0.001 


-0.001 



k 7 =-7t 



outer e-FS 



(d)inner e-FS (FS5) 
k 7 =;i 



k 7 =-7i 




FIG. 9: (Color online) Obtained SC gap functions for g = 
0.204 and U = 1.023. (a),(b) SC gap functions in the e-FSs 
on k z — and k z — n planes. (c),(d) 3D gap functions on 
the outer and inner e-FSs. The green lines represent the gap 
nodes. 

As increasing the value of U (or reducing ni mp ) fur- 
ther, the sign of the SC gap on the outer e-FS is com- 
pletely reversed, and small closed loop-nodes appear only 
on the inner e-FS. The obtained SC gaps are nearly s±- 
wave state. Figures 9 (a) and (b) show the obtained gap 
functions on the e-FSs for g = 0.204 and U = 1.023 
(a c = 0.980, a s = 0.869). They are approximately given 
by shifting the gaps in Figs. 8 (a) and (b) downwards. 
The obtained eigenvalue is Ae = 0.50. Figures 9 (c) and 
(d) show the obtained 3D gap functions on the outer 
and inner e-FSs, respectively. Apparently, closed nodal 
loops appear in the high curvature part on the inner e-FS, 
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whereas no nodes appear on the outer e-FS. This numer- 
ical result is consistent with the recent ARPES measure- 
ment by Yoshida et al [33]. On the other hand, the SC 
gaps on the h-FSs are similar to those in Figs. 7 (a)-(c). 

In Figs. 7-9, we fixed the impurity parameters as 
n-imp — 0-03 and 1=1. Now, we discuss the SC gap 
functions for general impurity parameters. Figure 10 (a) 



shows the U- 



phase diagram for both 1=1 and 



I = 0.3. The solid (dashed) lines represent the bound- 
aries between s++ wave and nodal-s wave, or nodal-s 
wave and s± wave for 1=1 (I = 0.3). As decreasing a s 
or increasing a c , the following crossover would be real- 
ized: (i) full gap s±-wave — > (ii) nodal s-wave — > (iii) full 
gap s ++ -wave. When both U and g are fixed, the same 
crossover occurs when rii mp increases. The residual resis- 
tivity for 1=1 derived from the linear response theory 
is about 20 /if2cm per rii mp = 0.01. 

We note that, in the present numerical calculation us- 
ing 3D model, line-nodes can appear even if 7ij mp = 
as shown in Fig. 10. In contrast, in the previous calcu- 
lation using 2D model [20, 42], we could not obtain the 
line- nodes for nj mp = 0, since the SC state changes from 
the s_)-+-wave to s±-wave discontinuously as U increases. 




7=0.3 
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U 



1.1 
1.0 
0.9 



(b)/=l \a s =l 




5+ wave 










nodal s 


^ ++ wave 




a c =0.98 : 







0.02 

"imp 



A 2 C =3 
A 2 C =2 



0.04 



FIG. 10: (Color online) U-nimp phase diagram for a c = 0.98 
obtained for (a) / = 1 and 0.3 with A c = 1 and for (b) 1=1 
with A c = \/2 and \/3. 

In the present study based on the RPA, the s++ <H> s± 
crossover is realized in case of a s <C a c (= 0.98) for 
"■imp <0.1. One of the main reasons would be the factor 
3 in front of V s in eq. (13), reflecting the SU(2) symmetry 
of the spin space. However, this factor 3 might be over- 
estimated since the recent polarized neutron scattering 
measurements indicates the relation xl(Q) ^ XL.yiQ) 
above T c due to the spin-orbit interaction XI ■ s [48, 49]. 



Moreover, we have recently improve the RPA by in- 
cluding the VC, and found that orbital fluctuations 
strongly develop in the Hubbard model [10]. Then, the 
orbital susceptibility is x c (q) = (X c (q) + x (q))(l — 
f c (X c {q) + x°(g))-\ where X c (q) is the charge VC for 
the irreducible susceptibility. According to Ref. [10], 
the magnitude of the three-point vertex is estimated as 
A c = 1 + X c (q)/x°(q) ~ 2, and then eq. (13) would be 
replaced with V c = iA^f c x c f c . Figure 10 (b) shows the 

C^-«imp phase diagram for A c = \[2 and -\/3, in case of 
1 = 1. We find that the s ++ -wave region is widely ex- 
tended, and the nodal-s-wave region is also widen. The 
obtained gap structure in the crossover regime for A c > 1 
is the loop-shape nodes shown in Figs. 7-9. 



IV. DISCUSSIONS 

In previous sections, we analyzed the gap equations 
based on the three-dimensional ten orbital model for 
BaFe2(As,P)2- When orbital fluctuations are solely de- 
veloped, fully-gapped s++-wave state is realized. On the 
other hand, when spin fluctuations are solely developed, 
we obtain the s±-wave state with horizontal node on a 
h-FS. During the crossover between s ++ -wave and s±- 
wave states due to the competition between orbital and 
spin fluctuations, the loop-shape nodes appear on the 
outer (inner) e-FS when the spin fluctuations are slightly 
weaker (stronger) than the orbital fluctuations. The ob- 
tained phase-diagram is shown in Fig. 10. We stress that 
all three h-FSs are fully-gapped during the crossover, 
since the SC gap on the z 2 -orbital originates from the 
inter-orbital nesting between different h-FSs. 

The crossover from the s±-wave state to the s++-wave 
state is also induced by increasing the impurity concen- 
tration. In this study, we considered the orbital-diagonal 
on-site impurity potential at Fe i-site, considering the 
Fe-site substitution by other elements. In this case, 
inter-band impurity scattering is always comparable to 
intra-band one, as shown by the T-matrix approxima- 
tion in Ref. [25]. For this reason, when the spin fluc- 
tuations are solely developed, the realized s±-wave state 
with T c o ~ 30K is suppressed by small amount of im- 
purities, with small residual resistivity po ~ 5z _1 fiflcm 
(z -1 = m*/m is the mass-enhancement factor). Since 
~ 3, we can safely expect that the SC state in dirty 
Fe-based superconductors (say po ~ 100/if2cm) would be 
the s ++ -wave state due to orbital fluctuations. 

There are many important future issues. As we dis- 
cussed in Sec. IIIC, one of our important future prob- 
lems is to study the present 3D ten-orbital Hubbard 
model (g = 0) using the SC-VC method developed in 
Ref. [10]. Using this method, we have recently shown 
that the s++-wave state is realized in the 2D five-orbital 
Hubbard model (g = 0). Note that the ferro-orbital 
fluctuations induced by the VC, which explain the or- 
thorhombic structure transition in Fig. 1, enlarge T c fur- 
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ther [52]. It was recently confirmed that the mechanism 
of orbital fluctuations due to VC had been realized even 
in a simple two-orbital model, using the SC-VC method 
as well as the recently developed two-dimensional renor- 
malization group analysis (RG+cRPA method) [43, 44]. 

Another important future issue is to include the self- 
energy due to orbital and spin fluctuations, E s and E c , 
into the gap equation. They are given as 



(24) 



where positive (negative) sign corresponds to £ = c (£ = 
s). The real and imaginary parts of the total self-energy 
= S c (fc) + S s (fc) represent the mass enhancement 
and quasiparticle inelastic scattering, respectively. Both 
effects suppress T c . Moreover, orbital and momentum 
dependence of would strongly modify the anisotropy 
of the SC gap functions. Thus, the self-energy correction 
in the gap equation will be important for the quantitative 
analysis of the SC gaps. 



V. SUMMARY 

In this paper, we studied SC gap structure using ten 
orbital model for BaFe2(As,P)2- When the orbital fluc- 
tuations due to inter-orbital quadrupole interaction (3) 
are strong, the s++-wave state is realized. In contrast, 
the s±-wave state is formed by strong spin fluctuations, 
mainly due to intra-orbital Coulomb interaction U. Both 
spin and orbital fluctuations would strongly develop in 
the optimally-doped regime near the O phase. In this 
case, we find that a smooth crossover between s++- and 
s±-wave states is realized by changing the interactions or 
impurity concentration, without large suppression in T c . 

During this s ++ f-> s± crossover, the loop-shape nodes 
are universally formed on the e-FSs, as a result of the 
competition between inter-orbital attractive interaction 
and intra-orbital repulsive interaction. This result is con- 
sistent with recent angle-resolved thermal conductivity 
measurement [34] and ARPES measurement [33]. Dur- 
ing the crossover, the SC gaps on the h-FSs are fully- 
gapped and almost orbital independent due to orbital 
fluctuations, consistently with recent ARPES measure- 
ments [32, 33]. 
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Appendix A: heavily electron-doped case 

In this paper, we have studied the gap functions for n = 
6.0. In this Appendix, we discuss the heavily electron 
doped case {n = 6.1). Figure 11 (a) shows the h-FSs on 
the k y — plane. As one can see, the xy-orbital h-FS 
disappears. 



(a) shape of the h-FSs (k y =0) 
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FIG. 11: (Color online) Obtained SC gap functions in the 
over-doped case (n — 6.1) for g — 0.25, (7 = and ni mp = 0. 
(a) fc z -dependence of the h-FSs, and (b) the obtained SC gap 
functions on the h-FSs. We also show the 3D gap functions 
on the (c) outer e-FS and (d) inner e-FS. 

Now, we consider the orbital-fluctuation-mediated 
s ++ -wave SC state. Here, we put g — 0.25 (a c = 0.98) 
and U = and obtained eigenvalue Ae is 0.87. First, we 
discuss the SC gaps on the h-FSs. Figure 11 (b) shows 
the k z dependence of the SC gaps on the h-FSs in the 
k y = plane. In contrast to the case of n = 6.0, the 
SC gap on the z 2 -orbital hole FS becomes smaller, while 
each SC gap on the h-FSs on any plane perpendicular to 
the fc z -axis is almost isotropic. Next, we focus on the SC 
gaps on the e-FSs. Figures 11 (c) and (d) show the 3D 
gap functions on the outer and inner e-FSs, respectively. 
As with the case of n = 6.0, the SC gaps on the e-FSs 
are nearly isotropic on any plane perpendicular to the 
fc z -axis. 

Next, we discuss the spin-fluctuation-mediated s±- 
wave SC state. Here, we put U = 1.37 (a s = 0.98) 
and g = 0. The obtained eigenvalue Ae is 1.68. Figures 
12 (a) and (b) show the gap functions on the h-FSs in 
the k z = and k z — 7r planes, respectively. In common 
with the results of n = 6.0, the SC gap function on the 
z 2 -orbital h-FS is almost zero since the spin fluctuations 
in the z 2 -orbital do not develop. Figures 12 (c) shows 
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FIG. 12: (Color online) Obtained SC gap functions in the 
over-doped case (n = 6.1) for g = 0, U = 1.37 and ni mp = 0. 
(a),(b) SC gap functions on the h-FSs in the k z — and 
k z — tt planes, (c) k z dependence of SC gap functions on the 
h-FSs for k y — 0. (d),(e) 3D gap functions on the outer e-FS 
and inner e-FS. 



the k z dependence on the SC gaps on the h-FSs in the 
k y = plane. It clearly shows the horizontal line node 
on FS3 near k z = n. 



Next, we discuss on the e-FSs. Figures 12 (d) and (e) 
show the 3D gap functions on the outer and inner e-FSs, 
respectively. Unlike the case of n = 6.0, line nodes appear 
on the e-FSs even when g — 0. This result is consistent 
with the analysis in Ref. [31]: nodal gap appears on the 
e-FSs when the xy-orbital h-FS disappears because of 
the competition of different spin fluctuations; Q = (tt, 0) 
on xz + yz-orbitals and Q — (tt,tt/2) on ccy-orbital. In 
the presence of xy-orbital h-FS for n w 6.0, on the other 
hand, spin fluctuations develop at Q = (jr, 0) in all d- 
orbitals. Then, fully-gapped s-t-wave state is realized. In 
this case, the competition of orbital and spin fluctuations 
induce the loop-shape nodes as discussed in Sec. III-C. 
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